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Abstract 

We investigate cosmological constraints on primordial isocurvature and tensor 
perturbations, using recent observations of the cosmic microwave background and 
the large scale structure. We find that present observations are consistent with 
purely adiabatic initial conditions for the structure formation under any priors on 
correlations of isocurvature modes, and upper limits on the contribution of isocurva- 
ture and tensor perturbations are presented. We also apply the obtained constraints 
to some specific theoretical models, axion isocurvature perturbation models and cur- 
vaton models, and give some implications for theoretical models. 



1 Introduction 



Recent cosmological observations, such as the cosmic microwave background (CMB) and 
the large scale structure (LSS) provide us information on primordial perturbations which 
seed the structure of the present universe. All observations suggest that the primordial 
fluctuation is almost adiabatic and scale-invariant [HIS]. Inflation is the most promising 
mechanism to generate the scale-invariant adiabatic fluctuation in the early universe. On 
the other hand, primordial isocurvature perturbations are also generally generated, along 
with the tensor perturbations, in the inflation universe. Many possible sources and mech- 
anisms generating isocurvature perturbations are known such as axion, curvaton scenarios 
[3|, H] and multi-fleld inflation models. 

Therefore it is expected that constraints on primordial isocurvature and tensor per- 
turbations give us some useful information to build realistic inflation models and models 
in particle physics. Thus, constraints on primordial isocurvature perturbations have been 
investigated by many authors [HI El [71 El [9l HOl HH [12] (For recent constraints we refer 
to fl3[ [TU [l5]). However, there have been few investigations on cosmological models 
with both isocurvature and tensor perturbations. This is partly because in most inflation 
models tensor perturbations are expected to be small when (especially correlated) isocur- 
vature perturbations are generated. However still some models predicts generation of both 
isocurvature and tensor perturbations [ini [IZj. From phenomenological point of view, it 
is worth checking whether cosmological observations are consistent with purely adiabatic 
initial conditions even if we consider both isocurvature and tensor perturbations. 

In this paper we investigate constraints on cosmological models with both isocurvature 
and tensor perturbations in light of cosmological observations of CMB and LSS. We use 
data from two recent cosmological observations, CMB temperature and polarization power 
spectra from WMAP 3-year result (WMAP3) and galaxy power spectrum from SDSS 
data release 4 of luminous red galaxy sample (SDSS DR4 LRG). We only consider models 
with one isocurvature mode along with adiabatic and tensor modes, which are simple but 
suggestive for various models predicting generation of isocurvature and tensor modes. We 
investigate the isocurvature mode by using three different priors on correlation between 
isocurvature and adiabatic modes; 1) uncorrelated, 2) totally correlated and 3) generally 
correlated models. The reason why we investigate uncorrelated and totally correlated 
models separately is that there are some simple models predicting deflnite correlations. For 
examples, the axion isocurvature perturbation model produces uncorrelated isocurvature 
mode and the totally correlated one is predicted in curvaton scenarios. 

The structure of the paper is as follows. In Section [2] we briefly review the general 
initial perturbations of the structure formation which includes the isocurvature and tensor 
perturbations and we also gives the parametrization used to constrain the isocurvature 
perturbations there. In section [3] we give some examples of models with both isocurvature 
and tensor perturbations which are based on inflation scenarios. In section [H we show the 
method to obtain the constraints on the isocurvature and tensor perturbations from the 
combined set of cosmological observations. In section [5] we present constraints on isocur- 
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vature and tensor perturbations from CMB and LSS. In section [6] we apply the obtained 
constraints on the isocurvature perturbation to some specific models; axion isocurvature 
perturbation models and curvaton scenarios. Section [7] is dedicated to conclusions and 
discussions. 



2 Initial perturbations for the structure formation 

Scalar perturbations are generally decomposed into five modes [18]; adiabatic mode (AD), 
CDM isocurvature mode (CI), baryon isocurvature mode (BI), neutrino isocurvature den- 
sity mode (NID) and neutrino isocurvature velocity mode (NIV). In the framework of the 
linear perturbation theory, each mode evolves independently and observables in present 
universe such as the CMB angular power spectrum and the matter power spectrum are 
predicted by initial amplitude of each mode and their correlation. 

In this paper we adopt the definition of initial perturbations for the structure formation 
in [18]. We use X/(k) for representing the initial perturbation of each mode. 

C (for AD) 

'ScDM (for CI) 
X/(k) = { Sb (for BI) (2.1) 

(for NID) 

^V. (for NIV) 

The right hand side of Eq. (12.11) is evaluated at the beginning of the structure formation. 
C is the gauge invariant curvature perturbation and iScdm, and Si, are the entropy 
perturbations of CDM, baryon and neutrino, separately. = V,^ — is the relative 
velocity perturbation of neutrino to photon, f^, is the fraction of the neutrino species in 
the energy density of the radiations. For more detailed definition of each isocurvature 
mode, we refer to [TH] . 

When we investigate observational constraints on various isocurvature modes in sec- 
tion [5], we consider only CI, NID and NIV modes. This is because the contribution of 
CDM and baryon isocurvature perturbations are brought together into isocurvature per- 
turbations of matter 

c _ ^CDM e , c (o 9^ 

'-'m — '-'CDM + {^■^) 

Thus, the constraint on Sb is easily obtained from that on iScdm- 

The auto and cross power spectra Vab{k) of primordial perturbations can be written 

as 

P/j(A:)5kk' = ^{Xj{krXj{k')). (2.3) 

Here subscripts I and J represent the adiabatic (AD) and four isocurvature modes (CI, 
BI, NID and NIV). We assume power spectra can be approximated as power- law: 

ruik) =Arj -] , (2.4) 
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where ko is a pivot scale and we consistently take k — 0.05 Mpc~^ in the rest of this paper. 

Throughout this paper wc consider models with only one isocTirvature mode besides 
adiabatic and tensor modes. This simplification enables us to capture what models are 
plausible to generate the initial fluctuations in the universe including iso curvature and 
tensor modes and the resultant constraints on isocurvature and tensor perturbations can 
be applied to many theoretical models based on inflation scenarios and particle physics. 
Since we have known that the primordial perturbations mainly consist of adiabatic pertur- 
bations, it is convenient to normalize the amplitudes of power spectra by the amplitude of 
the auto power spectrum of the adiabatic mode, A\d. Thus, we parametrize initial power 
spectra for the scalar perturbations as 

A A ( ^ BaCosea\ s 
AiJ-A^r.[B^,^,e. Bl J' (2.5) 

where 



Ba = jAaa/AAD, (2.6) 



COsOa = AAD,a/^AaaAAB- (2.7) 

Subscripts a represent some isocurvature mode being considered (CI, NID or NIV). 

The tensor perturbations are also generated in inflation models. When we refer to 
tensor to scalar ratio r, we usually consider cases that the scalar perturbation is purely 
adiabatic. Since we are considering isocurvature perturbations along with the adiabatic 
perturbation here, we redefine r as 'tensor to adiabatic ratio'. The power spectrum of 
tensor perturbations Vg{k) is written as follows: 

n(^) = ^.f^j , (2.8) 

(2.9) 

where Ug is the spectral index of the tensor perturbation. 

Since each mode evolves independently within the framework of the linear perturbation 
theory, we can decompose perturbations of the fiuids by initial modes. As for CMB, the 
brightness function of photon is written as 

ei{k,n)^^Qi{Kv), (2.10) 
I 

where 0[ are brightness functions which evolve from different initial perturbation modes. 
We introduce a transfer function of photon for each mode T/(A;, 77): 

eiik,rj)^T/ik,ri)Xiik). (2.11) 
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Then we obtain the angular power spectra of CMB C;, 



Aad [C^'"' + cos 9aCr + BlCT + rCf 



(2.12) 
(2.13) 



where Q's are the angular power spectra in cases that the amplitudes of the initial pertur- 
bations Aij are set to be unity and subscripts adi, iso, cor and tens represent adiabatic 
auto, isocurvature auto, adiabatic-isocurvature cross and tensor power spectra, respec- 
tively. Q's are given by 



^iso 
^cor 



An 


r dk 


21 + 1 J 


T 


An 


r dk 


21 + 1 J 


T 


Att 


r dk 


21 + 1 J 


~k 



™ad,ad-1 



n.a,a-l 



"^AD.a — 1 



Tr{k)mk). 



As for the matter power spectrum, P{k) can be written in the same way, 
P{k) = Aad {P'^Kk) + 2B„cos^„P^"(A;) + BlP'^^k)] , 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



where the hatted power spectra P{k)''s are auto and cross power spectra with Au being 
unity. 



3 Isocurvature perturbation based on inflation 

So far we have considered generic models with isocurvature and tensor perturbations. 
In this section we consider a model with isocurvature and tensor perturbations based 
on inflation models. We assume there are two scalar perturbations generated during 
inflation. One is a curvature perturbation and the other is a isocurvature perturbation 
»S*. The curvature perturbation raises only adiabatic mode at the beginning of structure 
formation, whereas the isocurvature perturbation 5* can generally produce both adiabatic 
and isocurvature perturbations at the beginning of the structure formation. Therefore we 
can write 

Here in the left hand side C, and Sa are initial curvature and isocurvature perturbations 
for structure formation, separately. T's are transfer functions which represent how initial 
perturbations for structure formation are generated from perturbations during inflation. 
Since the curvature perturbation at over-horizon scale stays constant in the absence of 
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isocurvature perturbations, 7^,^, = 1. On the other hand 7^,5, and ^5^,5, depend on 
models. Then the initial power spectra for ( and Sa are given by 



Pad,ad(A;) = T^,c.ikfVcM + T^,sM^T^sM, (3-19) 

VADAk) = r^,sXk)rs^,sXk)VsXk). (3.20) 

VaAk)=%.,sSkfVsSk) (3.21) 



Generally C* and 5* may be correlated, such in the case of multi-field inflation models [161 
[T7] . However, we assume that perturbations C* and 5* are uncorrelated in this paper. 
Furthermore, we take power-law forms for the terms in the right hand sides of Eqs. (13. 19^ - 
(JMH) as 

^C,C.(^)'n.W=^adii(^^j , (3.22) 

r^^sMVsSk) = A^^2[j] , (3.23) 

(. \ ncoi — 1 

rs^,sXkfVsSk) = A,,,(^^^ , (3.25) 
where Acor and ricor are given by 

Acor = ±V'Aadi2^iso, (3.26) 

''^adi2 ~l" ''^iso /r, 

?^cor = ■ (3.27) 

Here the sign in the right hand side of the first line comes from a factor T^^s^Tsa,s* 
in Eq. (13.241) . which can be either positive or negative. Furthermore, we use following 
parametrizations: 

(3.28) 
(3.29) 
(3.30) 

Assuming single-field slow-roll inflation and iS* is the isocurvature perturbation for 
some scalar field (7^ inflaton) whose mass is negligibly light compared with the Hubble 
parameter during inflation, the following inflation consistency relations should be satisfied: 



^AD 


= ^adil + ^adi2 




= ^iso/^AD; 


Ba COS 6a 


= ^cor/^AD- 



We finally obtain the power spectra for CMB and matter, 

Ci = Aab [sin^ OaCj''^'^ + cos^ ^aCT'^'^ + Ba cos OaCr + BlC'r + rCi 
P{k) = Aad [sin^ ^^P^^dii + cos^ ^^Pf^'^ ^ e^P^°'{k) + BlP'''"{k) 



tens 



(3.32) 
(3.33) 



Note that there are two terms for adiabatic modes in each Ci and P{k). CMB and matter 
power spectra with subscript adil come from and those with subscript adi2 come from 
i5* . 



4 Analysis method 

We consider the fiat ACDM model, and take the standard value 3.04 for massless neu- 
trino species. We do not consider runnings in the spectral indices for scalar and tensor 
perturbations. 

Since we are considering isocurvature and tensor perturbations, there exist six extra 
parameters (P^, cos6'a, nAD,a, '^a.a, r, Ug) or (P„, cos6'a, nadi2, ^iso, Ug) that are 
absent for a purely adiabatic case. However, for obtaining sensible constraints from the 
present cosmological observations, it is not suitable to treat all these parameters as free 
parameters. In this paper we adopt some simplifications and fix the spectral indices 
{nAB,a, na,a,ng) to somc values. For the case considered in the previous section, we adopt 
the inflation consistency relation Eq. (13.311) to fix (nadi2 7 ''^iso, %)• These simplifications 
reduce extra parameters to three primary free parameters {Ba, cos 6a, r). 

Thus our models have the following nine primary parameters: 

(t^b, t^CDM, ^'sound, ^AD, ^adi, ^a, COS 6a,r). (4.34) 

We investigate three models separately depending on the correlation of isocurvature modes; 
1) uncorrelated, 2) totally correlated and 3) generally correlated isocurvature models. 
When we investigate uncorrelated (cos Oa = 0) and totally correlated isocurvature (cos 6a = 
±1) models, we fix cos^^ and when we investigate generally correlated isocurvature modes, 
we assign fiat prior probabilities on the cos^q in the range [—1, 1]. 

The likelihood of a model is assessed using the WMAP three-year (WMAP3) data 
and likelihood code [19], [20] and SDSS data release 4 luminous red galaxy sample (SDSS 
DR4 LRG) [2j. We include the nonlinear corrections for the matter power spectrum 
|21] . and analytically marginalize over a bias parameter b and a parameter for nonlinear 
correction Q^i. We modify the CAME code |i22j to generate CMB and matter power 
spectra. Likelihood surfaces are explored by Markov Chain Monte Carlo methods using 
CosmoMC [23]. We generate six chains for each model with isocurvature modes and their 
correlation and apply the Gelman and Rubin convergence test [21]. We finally obtain at 
least 150,000 samples for each model, and in some cases over 400,000 samples. 
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Figure 1: 1- dimensional likelihood distributions for uncorrelated isocurvature models with 
CI (left), NID (middle) and NIV (right) mode, respectively. In each panel we show the dis- 
tributions for models without tensor modes using WMAP3 data only (red full), with tensor 
modes using WMAP3 data only (green dashed), without tensor modes using WMAP3 data 
combined with SDSS DR4 LRG data (blue dotted), with tensor modes using WMAP3 data 
combined with SDSS DR4 LRG data (magenta dot-dashed). 



5 Constraints on isocurvature and tensor perturba- 
tions 

5.1 Constraints on the uncorrelated isocurvature models 

Firstly we investigate uncorrelated isocurvature models (cos 6'a = 0). For the uncorrelated 
isocurvature models with tensor mode we impose inflation consistency relation 

Ug = -r/8, (5.35) 

which is realized in a single-field slow-roll infiation models. 

We present Id- marginalized likelihood distributions for the uncorrelated isocurvature 
models in Figure [H We also show 95% confidence limits (c.l.) on Ba and r for each models 
with (without) tensor modes from the combination of WMAP3 and SDSS DR4 LRG data 
in Table [H The CDM isocurvature (CI) and neutrino isocurvature density (NID) modes 
are rather tightly constrained and there is no improvement in minimum x^- On the other 
hand, presence of neutrino isocurvature velocity (NIV) modes tends to be favored by the 
present CMB and LSS data, though not yet at decisive level. Thus, we find no statistical 
support for finite contribution from uncorrelated isocurvature modes and CMB and LSS 
data are consistent with purely adiabatic initial scalar perturbations. 

"^^By single-field inflation model we mean inflation model where the vacuum energy is determined by a 
single field and does not depend on other light fields. 
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CI NID NIV 




Figure 2: 68 % (full) and 95% (dashed) 2- dimensional constraints on the uncorrelated 
isocurvature models CI (left), NID (middle) and NIV (left) with tensor modes. We present 
constraints using WMAP3 data only (green) and combined with SDSS DR4 LRG data. 



We also find 95% limits on tensor modes. Comparing with the constraint r < 0.30 (95% 
c.l.) for the model with purely adiabatic scalar perturbations [2J, we find that the upper 
limits on r for models with uncorrelated isocurvature are roughly same as that for the 
purely adiabatic model. This is because uncorrelated the isocurvature modes (except for 
uncorrelated NIV mode) and tensor mode contribute to the large scale anisotropy of CMB 
positively and there are no parameter degeneracy. For uncorrelated NIV mode, situations 
are little different since CMB power spectrum for NIV mode is relatively similar to that 
for AD mode [ see, e.g.. Fig. 1. in |25]]. Thus, the upper limit on r for uncorrelated NIV 
models is higher than those for other isocurvature models but it is still comparable with 
that for the purely adiabatic model. 





CI 


NID 


NIV 


r < 


0.31(0.33) 
0.26 


0.51(0.54) 
0.25 


0.69(0.62) 
0.31 


Amin 


0(0) 


0(0) 


-l(-l) 



Table 1: Constraints on Ba and r at 95% c.l. for uncorrelated isocurvature models with 
tensor modes (without tensor modes) from WMAP3+SDSS DR4 LRG. We also show the 
changes of the minimum values from the purely adiabatic model. 

It is known there are some parameter degeneracies among fractions of isocurvature 
modes, Ba, and other cosmological parameters Uf,, ujqdm and riadi- These degeneracies 
are understood by recognizing that the constraints on isocurvature modes rely mainly on 
the angular scale and the peak hight of the first acoustic peak in the CMB TT power 
spectrum. The relative hight of the first acoustic peak to the anisotropy at large angular 
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Figure 3: 1- dimensional likelihood distributions for totally correlated isocurvature models . 
Considered isocurvature and tensor modes and combinations of data are same as Figured 



scale, increases as the baryon density increases through compressions of the photon-baryon 
fluid. It also increases as the CDM density decreases and the early integrated Sachs- Wolfe 
effect is enhanced. Finally the peak height increases as the spectral index increases which 
leads to larger primordial fluctuations in small scales. For CI and NID modes, increase in 
Ba decreases the relative peak hight of the acoustic peak. Thus some cancelations exist 
among Ba, ujb, ^cdm and and parameter degeneracies arise. But for NIV mode, increase 
in Ba does not decrease the peak hight much and parameter degeneracies in NIV models 
are weaker than in uncorrelated CI and NID models. Though the peak hight has also 
strong dependence of the optical depth r, the polarization power spactra (TE and EE) 
of WMAP3 constrains r tightly and no parameter degeneracy between Ba and r is seen. 
Some parameter degeneracies such as degeneracy between Ba and ujcdm are broken by 
inclusion of LSS data and the constraints on Ba improve (for CI and NID modes seen in 
Figure [T]) . 



5.2 Constraints on totally correlated isocurvature models 

Next we investigate totally correlated isocurvature models (cos6'a = ±1)- Firstly we 
define parameters for collecting both positively (cos6'a = 1) and negatively (cos^^„ = — 1) 
correlated isocurvature models as 

w u a / ^" (for cos6'a = 1) . , 

Ba = Ba COS ea = | (^^^ COS = -1) " ^^'^^^ 

B'a take either positive, and negative values. For totally correlated isocurvature models 
with tensor modes we assume scale invariant tensor perturbations Ug = 1. 

We show Id-marginalized likelihood distributions for totally correlated isocurvature 
models in Figure [31 95 % limits and relative changes in minimum from the purely 
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Figure 4: 68 % and 95% 2-dimensional constraints on the totally correlated isocurvature 
models with tensor modes. Contours are same as Figure El 





CI 


NID 


NIV 




< 0.056(0.087) 


< 0.118(0.173) 


< 0.130(0.101) 




> -0.129(-0.080) 


> -0.151(-0.090) 


> -0.189(-0.174) 


r < 


0.49 


0.44 


0.30 


A mm 


0(0) 


0(0) 


0(0) 



Table 2: Constraints for totally correlated models with tensor modes (without tensor 
modes) from WMAP3+SDSS DR4 LRG. We show B'^ and r at 95 % c.l. and the changes 
of the minimum values from the purely adiabatic model. 



adiabatic modes are presented in Table [21 We find no improvement in values and 
the present observations of CMB and LSS are consistent with the purely adiabatic initial 
conditions. For any isocurvature modes, the limits on Ba for totally correlated models 
are found to be more stringent than for uncorrelated models. This is because, for totally 
correlated models, the correlation terms Cf^^ and P^"^{k) give significant contributions 
to the CMB and matter power spectra (see Eqs. fl2.13p and fl2.17p ). which is not present 
for uncorrelated models. Therefore both CMB and matter power spectra are affected 
much if totally correlated isocurvature perturbations are present and limits becomes more 
stringent. 

On the other hand the constraints on tensor modes r for CI and NID modes becomes 
weaker, compared to those for uncorrelated models. This is because anti-correlated CI 
and NID modes decrease the anisotropics in large angular scales of the CMB TT power 
spectrum, which can be partly canceled by contributions from tensor modes. However, 
the constraints on tensor modes are not affected much for correlated NIV mode. 
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5.3 Constraints on generally correlated isocurvature models 




Figure 5: 1- dimensional likelihood distributions for generally correlated isocurvature mod- 
els with CI (left), NID (middle) and NIV (right) mode. In each panel, we show distributions 
for models without tensor modes (red full), with tensor modes imposed inflation consis- 
tency relations (green dotted) and with tensor modes with fixed spectral index Ug = 1 
(blue dashed). 



CI NID NIV 




Figure 6: 68% (solid) and 95% (dashed) 2-dimensional constraints on the generally corre- 
lated isocurvature models. We present constraints for models with CI (left), NID (middle) 
and NIV (right) mode, separately. In each panel, we show constraints on models without 
tensor modes (red), with tensor modes imposed inflation consistency relations (green) and 
with tensor modes with fixed spectral index Ug = 1 (blue). Black dashed lines represents 
cos = ±1. 



We finally investigate generally correlated isocurvature models (— 1 < cos 6a < 1). 
For models with generally correlated isocurvature modes and tensor modes we consider 
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Figure 7: 68 % and 95% 2-dimensional constraints on the totally correlated isocurvature 
models with tensor modes. Contours are same as Figure [61 



two versions of models, models on which the inflation consistency relations Eq. fl3.3ip is 



imposed |^ and models with flxed spectral index Ug = 1 for tensor modes. 

We present the Id-marginalized likelihood distributions for generally correlated isocur- 
vature models in Figure [5], and 2d likelihood contours in Figure [61 We also show 95% 
confldential limits for Ba and r, mean values and 68% confldential limits for cos^a and 
changes in minimum values from the purely adiabatic models in Table [SB Still, we 
flnd that observations are consistent with adiabatic initial conditions. 

The upper bounds for Ba are similar to the uncorrelated isocurvature models. These 
results are also guessed from the results obtained in Section 15.11 and 15.21 since, as we have 
seen, the allowed contribution of isocurvature perturbations are higher for uncorrelated 
isocurvature models than those for totally correlated modes. 

We also present constraints for isocurvature and tensor perturbations in Figure [71 We 
can see that when we impose inflation consistency relations, upper bounds for tensor modes 
r are roughly same as those for uncorrelated isocurvature models. This can be understood 
as follows. When correlations of isocurvature perturbations with adiabatic perturbations 
are either positively or negatively large (cos^6'a — 1), the spectral index for tensor modes 
Hg takes large negative values for fixed values for r, resulting in too much fiuctuations for 
CMB anisotropy at large angular scales, which is disfavored from observations. Therefore 
large r is allowed only when the correlation of isocurvature perturbations is small and 
hence the resulted bounds on r are similar to those for uncorrelated isocurvature models. 
On the other hand, when we take the fixed spectral index for tensor modes, = 1, 
correlations of isocurvature perturbations can become large and the bounds for r weaken 
as in the cases for totally correlated isocurvature and tensor perturbation models with 



'^^Since Eq. p.3ip assumes slow-roll inflation, we must care that samples in MCMC chains should not 
take large values for slow-roll parameter e = —ng/2 = — r/16(l — cos^ Oa)- However, we have checked that 
e for each sample takes no more than 0.1 and our use of Eq. (|3.3ip is consistent. This is because the CMB 
power spectrum at large angular scale disfavors such negative large value of Ug. 
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Hg = 1. 





CI 


NID 


NIV 


Ba < 


0.28 


0.31 


0.58 


COS 9a 


-0.25 ±0.52 


0.03 ±0.62 


0.05 ±0.50 


r < 


0.32 


0.28 


0.31 


Amin 








-2 



Table 3: Constraints for generally correlated models with tensor modes imposed inflation 
consistency relations on from WMAP3±SDSS DR4 LRG. We show 95% c.l. for Ba and 
r, mean values and 68% c.l. for cos6'a, and changes of the minimum values from the 
purely adiabatic model. 





CI 


NID 


NIV 


Ba< 


0.29 


0.41 


0.76 


COS 9a 


-0.04 ±0.43 


-0.10 ±0.47 


0.06 ±0.34 


r < 


0.43 


0.50 


0.73 


A mm 








-2 



Table 4: Same as Table [3] except for Ug = 1. 





CI 


NID 


NIV 


Ba < 


0.33 


0.47 


0.59 


COS 9a 


-0.06 ±0.34 


0.06 ±0.45 


0.14 ±0.40 


Ay2. 

Amm 








-2 



Table 5: Constraints for generally correlated models without tensor modes from 
WMAP3±SDSS DR4 LRG. 



6 Application 

In this section we apply the constraints on the isocurvature perturbation obtained in the 
previous section to some specific models of particle cosmology. We investigate two kinds 
of models, axion isocurvature perturbation models and curvaton scenarios. 
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6.1 Constraints on axion isocurvature perturbation and inflation 
models 



Axion, which is originally proposed as a remedy for strong CP problem in QCD [261 [27| [28] . 
is a candidate for CDM. The properties of axion, such as the decay constant and its 
couplings to ordinary matters are constrained from various observations of astrophysical 
and cosmological phenomena [29l [30]. In inflationary universe, the axion field has CDM 
isocurvature perturbations JSV, ^32), '331 EH |351 [36l [371 EH] and they are constrained from 
observations of CMB and LSS [39, .40j. Firstly, we briefly review how the axion becomes 
CDM and its isocurvature fluctuation arises in the early universe. 

We consider the case where the PQ symmetry is spontaneously broken when the uni- 
verse is at the stage of inflation. During inflation the expectation value of the axion field 
is very smooth but fluctuates by the amount of the Hubble parameter H^f. The mean 
value of the axion x and its fluctuation 6x can be represented as 

X = faOi, (6.37) 

5x = H,^i/27T. (6.38) 

Here, is the axion decay constant and 6i is the initial phase of the axion field which 
takes an arbitrary value between — vr and vr. 

When the cosmic temperature is much higher than the QCD scale (T ^ Aqcd), the 
axion has no potential and its field value stays constant. As the universe expands and its 
temperature decreases, the universe undergoes the QCD phase transition and the axion 
obtains mass which depends on the temperature T as [H] 




m^{T) = Xm^{T = 0)1- , (6.39) 



where A ~ 0.1 and p ~ —4. When the axion mass becomes equal to the Hubble parameter 
[m^(T) ~ H{T)) ], the axion field starts to oscillate. After the axion starts oscillation its 
energy density scales as and behaves as CDM. The density parameter of the axion is 
given by 

where 7 is the dilution factor. If there occurs no entropy release after the axion starts 
oscillation, 7 = 1. The mass of the axion at zero temperature is determined by its decay 
constant fa 02] as 

m,(T = 0) = 1.3xlO-w(^^A_j . (6.41) 
Thus Eq. fl6.40l) can be rewritten in terms of fa as 



1.0 X 10-^ X 7gf ( "^Q^" ) ^'U—^ V^'. (6.42) 

^ ' V200MeVy I lOiOGeVy ^ ^ 
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The axion isocurvature (entropy) perturbation is written as 

5n^ 6n. 



(6.43) 



where and are the number densities of axion and photon, respectively. Axion 
isocurvature perturbation is given by the fluctuation of the axion field during inflation, 

5, = 2^ = (6.44) 

where we have used Eqs. (16.371) and (16.381) at the second equality. 

We consider the general case where CDM consists of axion and other particles and 
assume that only axion contributes to isocurvature perturbation. Then the CDM isocur- 
vature perturbation are given by 



>CDM 



l^CDM 



S^. (6.45) 



The curvature perturbation ( and tensor perturbations /i+,x are also generated during 
inflation and is written as 



where is the field value of the inflaton and Mpi = y/SrcG is the reduced Planck mass. 
We thus obtain power spectra of the adiabatic, CDM isocurvature and tensor modes as 

= A 

^AD - 1 = -6e + 2r], (6.51) 
nci-l = ng = -2e. (6.52) 

Here we assume slow roll inflation and e and rj are slow roll parameters given by 

.^iMj,(^^)\ (6.53) 
, - Ml.^. (6.54) 
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Figure 8: Constraints on the axion decay constant and Hubble parameter in the inflation 
universe. The colored regions are excluded by cosmic density of the axion (red), SN1987A 
(blue) and axion isocurvature perturbation (green). 



Since the axion isocurvature and the curvature perturbations are uncorrelated, cos ^ci = 
0. Using Eq. with Eqs. ([635]), fICTD and flCTD . we obtain Bci and r as 

u;, 2v/2iMpi 

Bci - 7^ , (6.55) 

^CDM Jafi 

r = 16e. (6.56) 

We are now ready to study constraints on axion and inflation models. From now on, 
we take Aqcd = 200 MeV, 6i = 1 and assume no entropy release occurs after axion starts 
oscillation, i.e. 7 = 1. Firstly we obtain a upper bound on the axion decay constant fa from 
the requirement that the energy density of axion should not exceed the observed matter 
density in the present universe, c<j^ < u;cdm- Combined with the lower bound obtained 
from supernovae 1987a [43J the axion decay constant should be in following range. 

10^°GeV < /a < 4.1 X lO^^GeV. (6.57) 

With using the bound for the CDM isocurvature mode in Table [T] we obtain the limits 
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on the inflation parameters as 



Hinf < lO^GeV, (6.58) 
e < 10"l^ (6.59) 
-0.05 < 7] < 0.06. (6.60) 

We also present obtained bound in the fa-Hmf plane in Figure [HI From Eqs. fl6.48p . fl6.40p 
and fl6.55p the ratio of the CDM isocurvature perturbation to the adiabatic one Bqi is 
written as 

(6.61) 

Therefore the upper bound on the Hubble parameter H^^f during inflation becomes lower 
as fa takes larger value. 

The resultant constraints on fa and ifinf are comparable with those in [ID], where the 
authors used the constraints on isocurvature and tensor modes derived by considering 
models with either of them, not both. We have analyzed models with both isocurvature 
and tensor modes but constraints on ifinf have not improved much. This is because the 
axion model predicts much less tensor perturbation than isocurvature one since fa -C Mpi. 
Thus, the obtained constraints on and i^inf do not change by inclusion of tensor modes. 
We can say oppositely that if nonzero contributions of tensor modes are suggested by 
future observations, the axion isocurvature model will be completely excluded. 

In the case where the initial misalignment of axion field is accidentally much smaller 
than its natural value, 6'j -C 1, the constraints on the axion decay constant weakens since 
initial amplitude for the oscillation of axion field becomes smaller. 



f'J'Ol (6.62) 



The constraints on ifinf also weakens. Although the amplitude of the isocurvature per- 
turbation in axion field becomes larger by decrease of initial misalignment, however, its 
fraction in CDM isocurvature perturbation becomes smaller since the fractions of axion in 
CDM becomes lower. 

Bci oc f'J'e,H,^f. (6.63) 

As we stated in the early part of this section, we have so far considered the case where 
the PQ symmetry is spontaneously broken during inflation. When the PQ symmetry is 
not broken during inflation or restored by the reheating after inflation, the inflation scale 
i^inf is not bounded by the constraints on the CDM isocurvature perturbations. 

6.2 Constraints on curvaton models 

In curvaton scenarios curvature perturbations are generated from the fluctuation of a 
scalar field (= curvaton) which is isocurvature at the epoch of inflation. We firstly briefly 
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review curvaton scenarios and then apply the constraints obtained in the previous section 
to them. 

We represent a curvaton field as a and an infiaton field as 0. Here we consider the case 
that the curvaton field is sufficiently light compared with the Hubble parameter during 
inffation. Then the mean value and ffuctuation of the curvaton field are given as 

a = ai, (6.64) 

5a = ^ (6.65) 

We represents the curvature perturbation generated during infiation as C* and the isocur- 
vature perturbation of the curvaton field as 5o- = 26a /(Xi. 

Until the Hubble parameter of the universe becomes below the mass of the curvaton 
mass, the expectation value of the curvaton field is constant. After the Hubble parameter 
becomes comparable to the mass of the curvaton, the curvaton field starts oscillation and 
its energy dominates the universe. When the curvaton starts dominating the density of the 
universe, its fiuctuation generates the curvature perturbation. After the curvaton decays, 
its energy turns into the radiation. If the curvaton produce the CDM, baryon or lepton 
number, their fiuctuations also obey the fiuctuation of the curvaton before its decay. Then 
various perturbations that are relevant for the structure formation are given by 

C = C + l^^, (6.66) 

ScBM = {rcBM - 1)5^ + fiy<S,y, (6.67) 

S, = (rB - l)S^ + US,, (6.68) 

S. = f(^^^ (rL-l)5., (6.69) 

where tcdm; i^b, are the fractions of the CDM, baryon number and lepton number 
densities produced by or after the decay of the curvaton in the present densities. ^ is 
the neutrino asymmetry parameter and we keep only the leading term of order in ^/vr in 
Eq. (16.691) since ^ is constrained from Big Bang Nucleosynthesis (BBN) [H] using the the 
observed helium abundance in HSl as 



1^1 < 0.07. (6.70) 

The neutrino isocurvature density perturbation in Eq. fl6.69p comes from the isocur- 
vature perturbation in lepton number density 5l = {5ni^/ni^—5n^/n^) [47J. This is because 
nonzero lepton number density in the universe nx, 7^ affects the energy density of neutrino 
via changing the distribution function function of neutrino through nonzero chemical po- 
tential. The lepton number density and neutrino energy density are both written in terms 
of neutrino asymmetry parameter ^ as 

31 



C(3)^3 




(6.71) 
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120' 




(6.72) 



Keeping only leading terms in ^/vr, we can relate the neutrino isocurvature perturbation 
Sy and the isocurvature perturbation for lepton number density iSl E1 

2 




(6.73) 



which yields Eq. (16.691) . 

More generally, the curvaton possibly decays before it completely dominates the uni- 
verse. We therefore phenomenologically parametrize the various perturbations by using 
= Pal Pti the ratio of curvaton energy density just before its decay to the total en- 
ergy density just after the curvaton decay, and then initial perturbations for the structure 
formation are written as 



#4 



C = C + 

^CDM = {f CDM — r)'5o- + fuSy 

Sb = (^B - r^)Sa + fuS^ 



(6.74) 

(6.75) 
(6.76) 

(6.77) 



In the case the curvaton decays after it completely dominates the universe and its energy 
turns into radiation nearly completely, = 1. 

Now we are prepared to obtain constraints on curvaton scenarios. For simplicity, 
we consider the case that the curvature perturbation generated at the inflation epoch is 
negligible (C* = 0) and the curvature perturbation is created by the curvaton. In that case 
the isocurvature perturbation is completely correlated with the curvature perturbation. 
Then, we can represent the parameters B'^ in Eq. (15.361) as 



B' 



'"cdmN 



CI 



B' 



BI 



B' 



NID 




28(1-/,) \ 



1-^ 



(6.78) 
(6.79) 

(6.80) 



"^^We simply assume there is no difference between perturbations in the temperatures of neutrino and 
photon. This is because photon and neutrino are thought to be coupled in the early universe at temperature 
0(1) MeV and their temperature keep fluctuating in the same way after the neutrino decoupling. 

"^^^ Authors in |48j used different parametrizations. Our parametrizations tr corresponds to Ar in [48j 
with taking Am = = 1. 
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Using the constraints on totally correlated isocurvature models without tensor modes 
obtained in Section 15.21 we obtain the following limits on BW^: 



- 0.029 < 1 - ™ < 0.027, (6.81) 
-0.134 < 1 - ^ < 0.133, (6.82) 

-7.2 X 10"^ < (5)^ (l - ^) < 3.7 X 10^^ (6.83) 

These constraints on curvaton scenarios are slightly stringent compared to those in [49] . 
and roughly same as those in [12] and [T3] . 

For = 1, the constraints imply that both CDM and baryon number should be created 
by or after the decay of the curvaton (rcoM — — 1). On the other hand, production 
of the lepton number is not constrained since no observation at present indicates the 
presence of non zero lepton number in the universe and ^ is consistent to zero. If we take 
a natural assumption that the lepton number should be comparable to the baryon number 
tzl/s ~ ub/s ~ 10~^° then the neutrino asymmetry parameter ^ should be of order 10~^. 
With such a small value of ^ the constraint Eq. fl6.83p then leads to 

< — < 10^ (6.84) 

rR 

and unless tr < 10~^ no restriction is assigned in generation of lepton number. Conversely, 
if nonzero fraction of the neutrino isocurvature density fluctuation is favored by future 
observations, the existence of large lepton number asymmetry may be suggested. 

We finally make a comment on the case where both and contribute to the initial 
perturbations for the structure formation. In this case, the constraints are weakened by a 
factor ~ >S<^/(3C*/rR + S^). 



7 Conclusion 

We have presented constraints on isocurvature and tensor perturbations from the combi- 
nation of CMB and LSS data. We have considered models with one isocurvature mode 
(CI, NID or NIV) and tensor modes. As for correlation of the isocurvature mode to the 
adiabatic mode, we have investigated three models; uncorrelated, totally correlated and 
generally correlated isocurvature models. 

For totally correlated isocurvature models, the contribution of isocurvature perturba- 
tion is severely limited Ba < 0.1 ~ 0.2. For uncorrelated and generally correlated isocurva- 
ture models we obtain Ba < 0.3 ~ 0.7 and upper limits are a few times larger than those 

"^^We here used the standard value for the neutrino fraction in the energy density of the radiation, 
fi, = 0.40. However if the large lepton asymmetry exists the thermal history of the neutrino is modified 
so that fi, is changed and the structure formation is also affected. We refer readers to [55] for various 
effects of the lepton asymmetry on the structure formation. Here we assume that the lepton number, if 
any, is sufficiently small and the thermal history of the neutrino is not affected. 
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for totally correlated models. Compared with other recent constraints on isocurvature 
models without tensor modes, our limits are roughly same even if contribution of tensor 
modes is included. 

We have also obtained the upper limits on the tensor mode taking the isocurvature 
mode into account. The limits are strongly depends on the isocurvature modes and its 
correlation included in the models. For CI and NID modes, constraints for uncorre- 
lated isocurvature models are similar to those for purely adiabatic models, but constraints 
weaken when correlation with adiabatic modes are included. For NIV modes, constraints 
for both uncorrelated and totally correlated models are similar as those for purely adiabatic 
models, but for generally correlated model, constraints loosen significantly. 

Finally we have found no significant improvement of fo^' models with isocurvature 
and tensor mode. Thus we conclude the initial conditions of the structure formation are 
still consistent with completely adiabatic ones. 

We have also applied the obtained constraints to some specific models which leads 
to the isocurvature perturbations, the axion isocurvature perturbation model and the 
curvaton scenario. Since the axion decay constant is bounded around 10^^ GeV, the scale 
of inflation i^inf which determines the amplitude of the axion fluctuation is constrained to 
be below 10^ GeV. Thus, very low scale inflation is required. As for the curvaton scenario, 
when the curvaton dominated the universe before its decay, we have shown that CDM and 
baryon number observed in the present universe should be created by or after the decay 
of the curvaton, otherwise too large isocurvature fluctuation is produced. However the 
generation of lepton number is not constrained by current cosmological observations. 
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